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ABSTRACT 

The Owens Valley mm- Array has been used to map the CO 1-0 emis- 
sion in the inner 2'. 5 of the grand design spiral galaxy M51 at 2" — 3" 
resolution. These new images reveal the molecular spiral arms with un- 
precedented clarity - the emission in the two major arms (NE and SW) 
originates from supermassive cloud complexes, Giant Molecular Asso- 
ciations (GMAs), which are for the first time resolved both along and 
perpendicular to the arms. The overall morphology of the CO emission 
is symmetric in reflection about the nucleus with major complexes oc- 
curing opposite each other in the two major arms. On the other hand, 
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the CO line flux in the area of the SW arm closest to the nucleus is ap- 
proximately twice as bright as that from the analogous location in the 
NE arm. 

Streaming motions can be studied in detail and appear with great 
clarity along the major and minor axes of M51. The streaming veloci- 
ties are very large, 60-150 km s _1 . Our maps offer, for the first time, 
sufficient resolution to resolve the structure in the molecular streaming 
motions. Both the radial and tangential velocity components show steep 
gradients in qualitative accordance with predictions by the density wave 
models of Roberts and Stewart (1987). Our data thus support the pres- 
ence of galactic shocks in the arms of M51. In general, velocity gradients 
across arms are higher by a factor of 2-10 than previously found. They 
vary in steepness along the spiral arms, becoming particularly steep in 
between GMAs. The steep gradients cause conditions of strong reverse 
shear in several regions in the arms, and thus the notion that shear is 
generally reduced by streaming motions in spiral arms will have to be 
modified. Of the three GMAs studied on the SW arm, only one shows 
reduced shear. 

We find an unusual structure, an expansion in the NE molecular arm 
at 25" radius (1.2 kpc) southeast of the center. This broadening occurs 
right after the end of the NE arm at the Inner Lindblad Resonance (ILR). 
Multiple-peak spectra, velocity twists and structures with apparently 
high velocity dispersion are associated with this feature. 

Bifurcations in the molecular spiral arm structure, at a radius of « 
3.3 kpc (73"), may be evidence of a secondary compression of the gas 
caused by, and occuring near, the 4/1 ultraharmonic resonance. 

Several molecular spurs protrude from the main spiral arms, in par- 
ticular on the western side of the map. Inside the radius of the ILR, 
we detect narrow (« 5") molecular spiral arms possibly related to the 
K-band arms found in the same region. We find evidence of non-circular 
motions in the inner 20" of M51. The magnitude of these deviations is 
20-30 km s _1 , and they are consistent with gas on elliptical orbits in a 
bar. 
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1. Introduction 

The Whirlpool Galaxy M51 at 9.6 Mpc 
(Sandage and Tammann 1975) is one of the 
closest grand design spirals, rich in molecu- 
lar gas, much of which is found in the spec- 
tacular spiral arms (e.g. Young and Scoville 
1983, Rydbeck et al. 1985, Lord and Young 
1990). The striking spiral pattern of M51 has 
stimulated many observational studies. Sin- 
gle dish CO studies of the molecular cloud 
distribution have been undertaken by Young 
and Scoville (1983), Rydbeck et al. (1985), 
Lord and Young (1990), Garcia-Burillo et al. 
(1993ab) and Nakai et al. (1994). Higher 
resolution aperture synthesis studies of the 
molecular spiral arms have been done by Vo- 
gel, Kulkarni and Scoville (1988), Rand and 
Kulkarni (1990), Tosaki et al. (1991), Adler 
et al. (1992) and Rand (1993a). Strong 
kinematic signatures of density waves are ve- 
locity discontinuities, or streaming motions, 
across the spiral arms (e.g. Roberts and 
Stewart 1987). The unusually large stream- 
ing motions (60-150 km s _1 ) of M51 im- 
ply a very strong density wave. The high 
resolution studies resolve massive cloud com- 
plexes, termed Giant Molecular Associations 
(GMAs) by Vogel, Kulkarni and Scoville (1988) 
along the main spiral arms and have also 
found a few apparently isolated complexes in 
the interarm region (e.g. Rand and Kulkarni 
1990; Adler et al. 1992; Rand 1993a). Possi- 
ble differences in properties between arm and 
inter-arm GMAs have been discussed (Rand 
1993ab) with the conclusion that the arm- 
GMAs are probably gravitationally bound ob- 
jects while the interarm GMAs are proba- 
bly unbound. It is still unclear, however, 
how density waves may affect the stability 
of the gas in the arms in detail. High res- 
olution studies of streaming motions within, 



and in between, GMAs have been lacking and 
it has been assumed that streaming does not 
change along the arm. The idea that arms 
are relatively shear-free environments where 
the streaming motions help cloud formation 
has been widely accepted. 

Spiral arms are more readily studied in 
gas, compared to stars (e.g. Lubow et al. 
1986; Lubow 1993). Sensitive CO observa- 
tions may therefore lead to a deeper under- 
standing of the underlying mechanisms be- 
hind spiral structure. Determining the loca- 
tion of resonances is particularly important 
since it helps distinguish between various the- 
ories of spiral structure. For instance, the- 
ories differ on whether spiral structure can 
extend beyond the 4/1 ultraharmonic reso- 
nance on the outside (e.g. Contopoulos and 
Grosb0l 1986; Artymowicz and Lubow 1992), 
and whether it terminates before reaching the 
Inner Lindblad Resonance (ILR) on the inside 
(e.g. Lubow 1993). 

In this first paper we will focus the discus- 
sion on streaming motions and spiral struc- 
ture. Newly discovered structure in the stream- 
ing motions on the main arms and the impli- 
cations for GMA formation and survival are 
discussed in section (4.1-4.3). The location 
of the 4/1 resonance and its possible relation 
to structures in the interarm emission is dis- 
cussed in section (4.4). Evidence for a molec- 
ular spiral in the inner kpc is presented in 
section (4.5) together with gas on non-circular 
orbits as indications of a central bar. Adopted 
properties of M51 are listed in table 1. 

In paper II we will focus on the smaller 
scale structure of the molecular gas. In par- 
ticular, the boundedness of clouds and possi- 
ble differences between on-arm, interarm and 
central region gas will be investigated, also in 
relation to tracers like Ha and infrared. 
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2. Observations 

Aperture synthesis CO mapping of M51 
was carried out with the Owens Valley Ra- 
dio Observatory (OVRO) millimeter array be- 
tween October 1995 and May 1996. SIS re- 
ceivers of the six 10.4m telecopes of the array 
provided typical system temperatures of 500 
K. The quasars 1150+497 and 0917+449 were 
used for phase calibration, and Uranus and 
Neptune for absolute flux calibration. Syn- 
thesized beams are 2."9 x 2."1 (PA=-81°) for 
robust (Briggs (1995)) weighting, and 3. "9 x 

3. "3 (PA=-43°) for natural weighting. 

Natural weighting, with all points cov- 
ered in the uv-plane treated alike, emphasizes 
shorter spacings, since the uv-tracks are more 
densely packed close to the origin. It gives 
the best signal-to-noise ratio, but at the cost 
of poorer beam shape. Uniform weighting as- 
signs identical weight to given cells in the uv- 
plane, thus maximizing the resolution at the 
cost of sensitivity (e.g. Sramek and Schwab, 
1989). Because of the limited uv-coverage 
common to mm interferometry, especially in 
mosaicing mode, this weighting scheme often 
yields poor results. Robust weighting, which 
we applied instead to produce our most highly 
resolved images, improves the signal-to-noise 
ratio over uniform weighting while still giv- 
ing higher resolution than naturally weighted 
maps. 

At wavelength A = 2.6 mm and resolution 
2. "5, a brightness temperature of ATb = 1 
K corresponds to 67 mJy beam^ 1 . The digi- 
tal correlator was centered at Ulsr = 460 km 
s _1 , and provided a total velocity coverage 
of 333 km s" 1 , from 351 km s _1 to 601 km 
s _1 . Data were binned at a resolution of 3 
MHz, corresponding to a spectral resolution 
of 7.8 km s" 1 . We observed a total of 19 



positions (listed in table 2) in three different 
array configurations (compact, low resolution, 
high resolution) — the shortest baseline was 
15m and the longest 240m. The largest struc- 
ture we could therefore image was 18". The 
primary beam is 65" (FWHM) and the grid 
spacing of the 19 positions was 30", chosen to 
give approximately constant noise across the 
map. 

We took 10 12-hr tracks per configuration 
for a total of 360 hours of telescope time. 
The observing scheme was designed to mini- 
mize calibration inconsistencies within the fi- 
nal map: All positions were observed sequen- 
tially, with only 2 minutes spent on each field 
at a time. Thus, each of the 19 positions was 
observed 10-15 times during each track, giv- 
ing a very homogeneous final data set. Indi- 
vidual maps were made and deconvolved with 
the NRAO AIPS and the Caltech DIFMAP 
(Shepherd et al. 1995) software packages. 
The AIPS LTESS task, a mosaicing procedure 
that corrects for primary beam attenuation, 
was used to combine all 19 positions into a 
single datacube. The rms noise is 19 mJy per 
channel (0.2 K) throughout the whole natu- 
rally weighted map (apart from the outer ~ 
3" before cutoff, where the noise reaches 35 
mJy per channel). 

3. Results 

3.1. Naturally weighted intensity map 

Figure la shows the naturally weighted 
integrated intensity map with 3. "5 resolu- 
tion. The molecular spiral arms appear with 
great clarity, and previously unknown struc- 
tures on, and between, the main arms are re- 
vealed. A comparison with optical pictures 
shows that the main CO arms trace the dust 
lanes beautifully. The map is dominated by 
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two main spiral arms of which the SW (Ml) 
arm is the brightest, its flux being greater by 
a factor of two than that of the NE arm (M2). 
The Ml arm is also broader than the M2 
arm by ^50%. The CO emission originates 
mostly from large cloud complexes (Giant 
Molecular Associations, GMAs). The bright- 
est (T B (CO) = 8 K) GMA on the Ml arm has 
a FWHM size of 11" by 6" (512 by 280 pc) and 
sits on a ridge of CO emission 20" x 7" in size. 
Along the SW Ml arm we count 16 GMAs, 
with sizes varying from 6" x 3.5" to 14" x 8". 
Local, distinct peaks in the intensity distribu- 
tion were counted as GMAs. Sizes for these 
peaks were determined using the AIPS IMFIT 
routine. Extended emission associated with a 
peak (at a level 20% or less of the peak) was 
not counted. The real number of GMAs in 
Ml is thus likely larger than 16. 

A more thorough analysis of the size and 
structure of the cloud associations will be 
given in paper II. Clouds previously studied 
by Rand (1993) (A4, A8 and A9) belong to 
the SW Ml arm and are indicated in figure 
la. The first gaps in intensity along the main 
arms occur symmetrically for both at a radius 
of 55" where, as figure la shows, the arms 
first cross the major axis. The NE spiral arm 
M2 can be traced inwards to a galactocen- 
tric radius of 28" (1.3 kpc), where the CO 
contours and the associated optical dust lane 
make a sharp turn eastward and the chain of 
GMAs ends. The first bright GMA on the 
SW Ml arm appears at a radius of about 22" 
(1.0 kpc) where its lowest contours bend east- 
ward towards the nucleus. At lower radii the 
cloud associations become poorly defined and 
fainter. Feature C (in figure la) is a somewhat 
peculiar structure, a broad expansion perpen- 
dicular to the molecular arm southeast of the 
center. There is a similar, but less spectac- 



ular, feature on the western side of the nu- 
cleus. Comparison of an optical image of M51 
with our CO map shows that the pattern of 
dust lanes spreads in a similar fashion to the 
CO emission, and that inside the C region the 
dust lane structure is disordered. Feature C 
occurs right after the ending of the main M2 
arm at the claimed position of the ILR (Tully 
1974b). We detect hints of molecular spurs, 
the features SI, S2, S3, S4, S5 in figure la, 
protruding at an angle (45° - 90°) from the 
main spiral arms mostly on the western side 
of the galaxy. 

A fairly recently discovered feature in M51, 
also apparent in our map, is the compact 
molecular gas concentration at the nucleus of 
the galaxy, first detected in HCN by Kuno et 
al. and mapped at a higher resolution of 1" 
in CO 2-1 by Scoville et al. (1997). 

There is evidence in the present study of 
bifurcations in the molecular spiral arm struc- 
ture - - possibly as a result of the 4/1 ul- 
traharmonic resonance caused by the spiral 
potential (e.g. Shu, Milione and Roberts 
(1973); Artymowicz and Lubow 1992 (AL)). 
Two faint chains of molecular clouds, Bl and 
B2 split off from the main spiral arms at a ra- 
dius of roughly 3.3 kpc (71"). Structures pre- 
viously classified as discrete interarm clouds 
(11,12 in Rand 1993a) are part of the Bl chain 
of clouds. 

3.1.1. The reality of the Bl feature 

The constant offset between the main arm 
Ml and the secondary arm Bl raises the con- 
cern that Bl is spurious, a sidelobe artifact of 
Ml, possibly the result of incomplete cleaning 
or miscalibration. S. Vogel (private commu- 
nication) raised this possibility because a new 
CO map of M51 made with the BIMA array 
does not show Bl, but when he simulated ob- 
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servations of the BIMA map using some (but 
not all of the OVRO configurations), an arm- 
like feature appeared that is similar to if not 
identical to Bl. In view of this possibility, we 
carefully reexamined evidence for and against 
the reality of the Bl arm. We conclude that 
it is probably real for the following reasons: 

i) The brighter emission peaks along Bl are 
clearly evident in IRAM and OSO single dish 
CO maps (e.g. Garcia Burillo et al. 1993a 
(fig 1), Rydbeck et al. 1991). The HI map 
by Rots et al. (1990) shows structure that 
is coincident with Bl. A number of faint Ha 
emission regions are also seen along Bl (Vogel 
et al. 1993). 

ii) A radial position-velocity cut through 
Ml and Bl reveals that although there is 
some overlap of the velocity ranges of emis- 
sion between Ml and Bl, the peak velocity of 
Bl is redshifted by 20 km s _1 relative to Ml. 
The velocity width of Bl is also narrower by 
a factor of two (see figure lc). It is significant 
too that, on the assumption of a flat rota- 
tion curve, Bl is at the expected velocity for 
its radius. Lastly, we note that the stream- 
ing velocity structure from the main arm to 
the bifurcated material is continuous at the 
beginning of Bl. 

iii) The new BIMA M51 CO map (Heifer 
et al. 1999, in preparation) does in fact show 
several low signal-to-noise features near the 
peaks in Bl. Since the BIMA map is of lower 
sensitivity and resolution than ours, we think 
that the two data sets may not be inconsis- 
tent. 

iv) While Vogel's simulation of the OVRO 
map does produce false arms, these do not 
agree very well in structure and position with 
our respective arms Bl and B2. The Bl 
feature in the simulation is shifted 5-10" in- 
ward from our Bl, and is straight, not curved. 



This last argument is of course not conclu- 
sive, since the simulation did not include the 
OVRO compact configuration (which proba- 
bly would reduce the spurious features) but 
in any case, if the residual sidelobes were the 
result of miscalibration, the simulation would 
never exactly reproduce the artifacts. 

v) The secondary structures survive deep 
cleaning, with various approaches such as 
varying gain and number of iterations, as well 
as the position and number of cleaning boxes. 
They are also visible in a robustly weighted 
map, where the sidelobes are less prominent. 
The B features are also apparent in a map not 
biased toward positive emission. 

vi) The Bl feature is well inside the pri- 
mary beam half power points (see table 2 for 
pointing positions) which means that it is not 
an edge effect of the map. The outer cut- 
off of the naturally weighted map (the outer 
dotted line in figure la) is somewhat (about 
5") inside the half power points of the outer 
pointings to show a map with (almost) con- 
stant noise. Note that the outer sensitivity 
cutoff of the robustly weighted map (figure 
lb) occurs at lower radius than for the natu- 
rally weighted one. 

Lastly, we note that the possibility that the 
B features are aliasing effects is excluded be- 
cause the size of the FFT was four times the 
map size. 

Between Bl and the bright Ml arm, there 
is a shorter chain of interarm clouds (feature 
D, indicated by a dotted line in Fig. la). 
These features also appear to be real struc- 
tures for the same reasons as above. 

3.1.2. Radial flux distribution 

The total integrated flux in the whole OVRO 
map is 2.7 x 10 3 Jy km s _1 . We have used the 
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OSO 20m single dish map of Rydbeck et al. 
(1985) to compare the total CO emission with 
that observed here interferometically. Inside a 
radius of 30" the flux recovered by the OVRO 
array is 810 Jy km s _1 while the OSO single 
dish finds 1600 Jy km s" 1 . Inside a radius of 
60" the OVRO flux is 1700 Jy km s" 1 and the 
corresponding value for OSO is 2800 Jy km 
s _1 , so that the OVRO array recovers 50-60 
% of the OSO single dish flux in the region 
studied. Comparing our map with the NRO 
single dish map of Kuno et al. (1995), we re- 
cover 30% of their single dish flux in the inner 
25". One large error source is thus the relative 
calibration among telescopes, with OSO and 
NRO obviously differing by almost a factor 
of two. Adler et al. (1992) claim to recover 
50% of the total OSO single dish flux in their 
BIMA aperture synthesis CO map. 

Figure 2a shows the brightness of each ra- 
dial bin, the OVRO mean radial CO distri- 
bution from the naturally weighted map. For 
the adopted, radially independent conversion 
factor (see table 1) between CO luminosity 
and H 2 mass, this figure shows the change 
in the mean gas surface density with radius. 
The gas surface density increases steeply be- 
yond a radius of roughly 8" and reaches a 
maxium (166 M pc~ 2 ) at a radius of 28" 
(1.3 kpc) which is that of the bright peak 
in the SW arm. There is then a sharp drop 
to a minimum at R — 52" (2.4 kpc), reflect- 
ing the first intensity gap in the main spiral 
arms. The brightness then increases some- 
what again, and after 90" it falls off steeply 
to the edge of the map at 95". From 28" to 
70", the average gas surface density drops by 
a factor of three. Kuno et al. (1995), from the 
NRO single dish map, find a drop by a factor 
of ^4 in the same interval, but their radial 
curve is less sharply peaked, either because 



of their poorer resolution, or because the flux 
undetected at OVRO is more smoothly dis- 
tributed. The single-dish gas surface density 
overall may drop somewhat faster with radius 
than that detected here, so it is possible that 
the interferometer recovers a larger fraction of 
the single-dish flux with increasing radius - 
consistent with an increasing arm-to-interarm 
ratio with radius. 

Figure 2b shows the variations of the CO 
intensity along the two major arms. The res- 
olution has been smeared to 10". The clumpy 
nature of the arms is evident, as well as there 
being a striking regularity of the distance, 
some 30" - 40", between the major concentra- 
tions of gas along the arms. These peaks even 
occur at (roughly) the same spiral phase on 
both arms. These symmetries were also noted 
by Garcia-Burillo et al. 1993a. A similar ef- 
fect in the optical was found by Elmegreen, 
Elmegreen and Seiden (1989) who interpreted 
it as an intrinsic part of the stellar wave pat- 
tern. 

3.2. Robustly weighted intensity map 

The robustly weighted map (figure lb) re- 
covers about 73% of the flux of the naturally 
weighted map in the same region, and much of 
the faint interarm structure is lost, although 
the brighter parts of the Bl and B2 features 
are present. The main arms are (not surpris- 
ingly) narrower in this map and the spurs are 
not as clearly visible. In the inner arcminute, 
less flux is lost (the map recovers about 85% 
of the naturally weighted flux) and structures 
that could be narrow spiral arms can be dis- 
cerned. 
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3.3. Comparison with previous aper- 
ture synthesis maps 

Because of the high sensitivity, high resolu- 
tion and improved uv coverage in our map, we 
find structures undetected by previous stud- 
ies. These structures include most of the 
spurs, the C feature, parts of the B and D 
features, and the narrow molecular arms in 
the inner arcminute. 

Our maps in general agree well with those 
obtained at OVRO by Rand and Kulkarni 
(1990, RK) at 8" resolution. There are hints 
in their data of the SI spur and the B2 bifur- 
cation; the positions of GMAs and emission 
gaps also agree with ours. There is some dis- 
agreement between the RK map and ours in 
the intensity distribution across the map; the 
SW arm is too faint in their map and Rand 
(1993, R93) also points out that that part of 
their map is suspect owing to calibration er- 
rors. R93 mapped two overlapping V fields 
that cover part of the Ml arm, and some in- 
terarm emission, with the OVRO array at 3. "5 
resolution; his arm structure and cloud posi- 
tions agree well with ours, although he misses 
much structure which we see (e.g. the A10 
cloud) owing to his poorer uv coverage. Rand 
found two interarm emission features (II and 
12 in R93) that we assign to the Bl cloud 
chain. Adler et al. (1992) presented obser- 
vations of M51 obtained with the Berkeley- 
Illinois-Maryland Association (BIMA) inter- 
ferometer; they observed two 2' fields with 
a synthesized beam of 7" x 11". Since their 
shortest baseline (8 m) was about half of ours, 
they were sensitive to larger structures. Given 
the much larger synthesised beam of the Adler 
et al. map, the qualitative agreement be- 
tween their map and ours is as good as can 
be expected. 



3.4. Velocity Field and Velocity Dis- 
persion 

The naturally weighted velocity field is 
shown in Figure 3 [Plate 1]. Velocities range 
from 350 (deep blue) to 550 (red) km s -1 , 
with the gas approaching us in the northern 
part of the galaxy and receding in the south. 
Streaming motions appear very clearly in the 
velocity field as jumps and discontinuities in 
the contours across spiral arms, in particular 
along the SW main arm, Ml. Radial stream- 
ing motions, best seen along the minor axis, 
are prominent in the Ml arm (west) while 
none are readily apparent in the east. Stream- 
ing motions along the major axis are more 
symmetrical with respect to the Ml and M2 
arms. 

Because of the significantly increased reso- 
lution of the robustly weighted map, stream- 
ing motions on the main arms are better re- 
solved than in the naturally weighted map, 
in particular where the Ml arm crosses the 
major axis. Sections of the robustly weighted 
velocity field are displayed in figures 6a and 
6b. 

The velocity contours associated with the 
C feature are somewhat twisted, indicating a 
complex velocity field in this region. Spurs 
SI and S4 show velocity reversals as they 
stretch out into the interarm region, S2 and 
S3 show strong velocity gradients but no re- 
versals, while S5 shows no velocity structure 
at all. Along S4 the velocity shifts from 500 
to 540 km s" 1 and then back to 500 km s _1 
as we move into the arm. Within SI the shift 
is from 400 km s _1 to 450 km s _1 and then 
back again. The red feature in the center of 
the map is the small rotating nuclear gas disk 
which, as previously noted by Scoville et al. 
(1997), is redshifted with respect to the sys- 
temic velocity by ~ 40 km s _1 . 
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The naturally weighted dispersion map is 
shown in figure 4 [Plate 2]. The map shows 
the one dimensional dispersion ranging from 
&id = 5 — 20 km s _1 . Only one value of aid 
is fit, even if the spectrum in question con- 
sists of multiple peaks. The red regions in 
the map indicate high dispersion, <7id^l8 km 
s _1 . Most of the high structures may be 
caused by unresolved streaming motions (e.g. 
where the western minor axis cuts the Ml and 
Bl arms). Positions of strong streaming mo- 
tions, indicated in the velocity field, do indeed 
correlate with regions of high dispersion. 

Other regions of disturbed velocity con- 
tours, such as the C region, and the bifur- 
cation areas also show spots of high disper- 
sion gas. In the C-region, several of these 
high dispersion spots are caused by multiple 
peaked spectra, but there is a small area of 
faint, broad line emission at the southern tip 
of the C feature, right before the feature nar- 
rows and <7id drops to 6-10 km s _1 . It is ev- 
ident too that the nucleus is a region of high 
velocity dispersion. 

Measured on the naturally weighted map, 
the average velocity dispersion in the Ml arm 
is < <7id >= 12 km s _1 and 7 km s _1 in the 
M2 arm before the emission gap at radius 55". 
After the gap, the average dispersion is about 
the same for both arms, < <7id >= 9 km s _1 . 
The average dispersion on the Bl structure is 
< Cid >— 6 km s -1 . The emission from the 
molecular spurs and the D clouds is in general 
of low dispersion (4-8 km s _1 ), except for a 
few positions along SI. The general result, 
that the dispersion of the interarm emission 
gas is lower than that of the arm emission, 
is consistent with the single dish result (at 
lower resolution) by Kuno and Nakai (1997) 
who suggest that the dispersion is higher on 
the arms because of the streaming motions. 



Indeed, < a ld > on the Ml arm is reduced 
from 12 to 9 km s _1 at the higher resolution 
offered by the robustly weighted map. 

3.5. Molecular mass 

We estimate the molecular mass in the in- 
ner 2/5 to be 1.9xl0 9 M , by adopting a stan- 
dard CO to H2 mass conversion factor (see ta- 
ble 1). Since between 70% and 40% of the to- 
tal flux is not detected by the interferometer, 
the real molecular mass in this region could be 
as high as 6 x 10 9 M assuming that the same 
mass conversion factor holds for the arm and 
interarm gas. However, according to Nakai 
and Kuno (1995), Arimoto, Sofue and Tsujo- 
moto (1996) and Adler et al. (1992), the CO 
to H2 conversion factor for the inner 4.5 kpc 
of M51 is lower by a factor of 2-3 than the 
Galactic mean value adopted here, increas- 
ing with radius. These values are based on 
CO-independent determinations of molecular 
masses by applying optical extinction mea- 
surements or assuming virial equilibrium and 
comparable to the lowest values suggested for 
our Galactic disk. 

3.6. Arm-interarm contrast 

Using the naturally weigthted map (figure 
la), we partitioned the observed emission into 
arm and interarm regions to determine the 
interferometric arm to interarm ratio. The 
arms were defined by the lowest contour in 
figure la. We have classified everything out- 
side the main stellar arms as interarm — for 
instance, the Bl and B2 features have been 
considered interarm as have the spurs pro- 
truding from the main arms. In Jy/beam km 
s _1 the mean flux in the arms is 3000 and 
that in the interarm regions is 190 resulting 
in an average arm-to-interarm ratio of about 
16. The interarm regions have a total inte- 
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grated flux of 390 Jy km s _1 , and the arm 
regions have 2300 Jy km s _1 . This gives a 
total flux ratio of 5.9, and with the standard 
conversion factor (see table 1) the mass in the 
interarm molecular gas observed here is about 
2.9 x 10 8 M . 

The actual arm-to-interarm flux ratio de- 
pends on how the missing flux is distributed. 
If the unseen emission is evenly distributed 
across the map, the actual arm-to-interarm 
ratio will be lower than that just derived. The 
total flux in our map is 2700 Jy km s^ 1 ; as- 
suming an equal amount is unobserved, and 
evenly distributed across the map, the arm- 
to-interarm contrast is reduced to about 3. 
The arm to interarm mass ratio would then 
be 1.3 instead of 6. In the same way, Vo- 
gel, Kulkarni and Scoville (1988) estimate 
from their interferometric CO map an arm-to- 
interarm mean flux ratio of 2.4-3.0. However, 
Adler et al. (1992) argue that when combin- 
ing a single dish map (the OSO 20m map of 
Rydbeck et al. 1985) with their aperture syn- 
thesis BIMA map, most of the missing flux is 
instead distributed in close association with 
the arms, making them broader, and they 
measure an arm-to-interarm mean flux ratio 
>4.6. 

Single-dish measurements detect all of the 
flux, but may lack sufficient resolution to 
clearly distinguish between arm and interarm 
regions. Previously measured single-dish CO 
arm-to-interarm ratios for M51 show average 
values of the mean arm-to-interarm flux ratio 
of 3-5, but also significant variation with ra- 
dius: Garcia-Burillo et al. (1993a) measure a 
peak arm-to-interarm ratio that increases with 
radius from 2 to 6. Kuno and Nakai (1997) 
find very similar values of 1.5-7, again increas- 
ing with radius, with the NRO telescope. 

The molecular arm-to-interarm flux ratio 



is interesting as an indication of the presence 
and strength of a density wave in a galaxy 
(e.g. Roberts and Stewart 1987; Elmegreen 
1988). The density contrast increases with 
increasing strength of the spiral gravitational 
field (e.g. Shu, Milione and Roberts (1973)). 
Roberts (1993), for example, finds an arm-to- 
interarm contrast of 6 for a 10% modal spi- 
ral driving. While our map is the first large 
scale interferometric data set to detect signif- 
icant interarm emission in M51, the contrast 
of 12 we find is still likely to be an upper 
limit. Thus, the arm/interarm contrast from 
our data cannot be used to constrain spiral 
density wave models beyond stating that it 
is generally compatible with them. Since the 
ISM passes into and out of the arms cyclically, 
the flow of matter must obey the continuity 
equation in order to conserve mass. Thus, if 
t arm is the time a molecule stays in the arm 
and ti n tcr is the time it spends in the inter- 
arm region, then ^ssn. — Mmtsi _ j n principle, 

farm ^intcr 

it should be possible to determine arm resi- 
dence timescales from this requirement. Our 
limits of arm-to-interarm mass ratio of 1.3 to 
5 are not restrictive enough to pursue this, 
but a combination with single dish data might 
make a detailed analysis feasible in the future. 

4. Discussion 

4.1. Streaming motions 

M51 is one of the galaxies where the strongest 
evidence for density waves has been found 
(e.g. Rand 1993b). The magnitude of the 
streaming motions is unusually high (e.g. Tully 
1974b; Rydbeck et al. 1985; RK; Vogel et 
al. 1993) amounting to about 100 km s _1 
in the plane of the galaxy - - a factor of 
two higher than predicted by linear theory, 
or non-linear simulations (e.g. Tully 1974b; 
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Roberts and Stewart 1987). In previous CO 
work at OVRO (Vogel, Kulkarni and Scov- 
ille 1988; RK; Rand 1993a (R93)) both radial 
and tangential streaming motions in the arms 
of M51 were detected, and streaming motions 
are also evident in single dish CO data (e.g. 
Rydbeck et al. 1985; Garcfa-Burillo et al. 
1993ab; Kuno et al. 1995). 

Density waves are expected to cause veloc- 
ity jumps and discontinuities across molecu- 
lar spiral arms. Tangential streaming is ex- 
pected to be the most apparent along the 
major axis, radial streaming along the minor 
axis. Along the major axis the highest local 
velocities (negative on the blueshifted side of 
the galaxy and positive on the redshifted side) 
are found on the outside of the arm, and on 
the inside one finds the lowest local veloci- 
ties. Outside of corotation, the situation is 
reversed. On the minor axis, a radial veloc- 
ity discontinuity, a decrease followed by an in- 
crease, across the arm is expected. In our map 
the density wave perturbations begin (for the 
Ml arm) at a radius of 22" (« 1 kpc), and it 
appears that the streaming motions are more 
distinct in the brighter, broader Ml arm than 
in the M2 arm. 

Roberts and Stewart (1987) (RS) have made 
simulations of the response of the interstel- 
lar medium to a spiral potential, treating in- 
dividual molecular clouds as point particles. 
In their model the spiral forcing amplitude 
relative to the central axisymmetric force is 
moderate, i.e. 10%. They studied mod- 
els that included dissipative cloud-cloud colli- 
sions and simplified cases without cloud-cloud 
collisions. They find that dissipative cloud- 
cloud collisions are necessary to stabilize the 
streaming motion pattern in time, and will 
produce stronger, sharper velocity and den- 
sity profiles than the non-collisional case. Fig- 



ure 5 shows the cloud density, tangential and 
radial velocity components plotted vs. spiral 
phase, from RS dissipative case simulations. 

Note that, for M51, the notion of a (quasi) 
steady density wave (e.g. Lin and Shu 1964) 
may be incorrect if the waves are triggered 
by swing amplification by the passing of the 
companion galaxy (e.g. Toomre 1981). In 
this case, the density wave pattern is transient 
and a steady streaming motion pattern may 
be unlikely. 

4.2. Tangential streaming 

According to the RS model (figure 5) the 
tangential streaming decreases as gas approaches 
an arm and then quickly increases across the 
arm. There is a very sudden change in slope, 
a sharply pointed trough, indicating a shock. 

Where both the Ml and M2 arms cross the 
major axis (at radius 55") there is a gap in 
CO intensity. We have therefore studied the 
streaming motions in these arms at ±10° off 
the major axis. Figure 6a and 6b show the ro- 
bustly weighted (2. "5 resolution) velocity field 
(contours) overlayed on the integrated inten- 
sity (grayscale) in the southern and northern 
arm-crossing regions, respectively. In figure 
6a the GMAs previously studied by R93, A8 
and A9, are marked as well as a GMA between 
these that we call A10. The non-linearity of 
the streaming motions is clearly evident here: 
The velocity contours are not evenly spaced 
across the arm, but instead a narrow region of 
densely spaced contours indicates a very steep 
increase in velocity. 

Furthermore, figure 6a clearly shows that 
streaming is highly uneven along the arm. 
Most of the streaming occurs before cloud 
A9 and through A10. In A8 some stream- 
ing can be detected but the gradient is not 
as steep as at A10. Streaming velocity gra- 
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dients through five GMAs and three inter- 
cloud regions (IC1,IC2,IC3 in figues 6a and 
6b) are listed in table 3. As is evident 
from the table, gradients vary significantly 
from GMA to GMA, the highest being that 
through A10, with 700 km s _1 kpc -1 and the 
lowest through A9 with 65 km s _1 kpc^ 1 . In 
general, gradients appear steeper in between 
GMAs. IC1 is the region upstream of A9 and 
between A9 and A10 and has the highest ve- 
locity gradients measured in this map, 800- 
1300 km s^ 1 kpc" 1 . Although Rand (R93) 
found a similar total magnitude of the stream- 
ing motions, he measured a significantly lower 
gradient of 100 km s _1 kpc -1 at the south- 
ern arm crossing. He also assumed the gradi- 
ent to be constant along the arm. The differ- 
ence probably results from Rand's poorer uv 
coverage and angular resolution. 

The highest total magnitude of the tangen- 
tial streaming motions, 150 km s _1 , is mea- 
sured through A10. The magnitude of the 
streaming motions is somewhat lower on the 
NE arm, about 100 km s^ 1 across A2. 

4-2.1. Arm shear 

Why do we see such structure in stream- 
ing motions and streaming motion gradients 
along the arm? Is this just an effect of the 
GMAs progressing across the arm, with A9 
perhaps having entered a region of low ve- 
locity gradient further downstream, or is it 
related to the cloud's internal response to 
streaming? Indeed, the presence of GMAs 
and the structure of the streaming are related, 
since the steepness of the gradient increases 
between the complexes. It seems that the 
presence of clouds does indicate lower shear 
regions, either because clouds form preferen- 
tially in such regions or because many of them 
have detached themselves more or less com- 



pletely from higher shear regions. Stream- 
ing motions have been suggested as conducive 
to molecular cloud formation. A local solid 
body rotation is supposedly created which 
increases shear time scales giving the gas 
enough time to form complexes before they 
pass into the harsh inter-arm region. The 
shear time-scale is the inverse of Oort's A- 
constant (e.g. Elmegreen 1988): 



where R is the radius and v c the circular 
velocity. It is clear from this relation that 
for the shear to be successfully counteracted, 
must be very close to ^, i.e. local con- 
ditions must approach solid body rotation. 
Rand (1993b) used this relation to estimate 
shear time scales. For clouds on the NE (M2) 
and SW (Ml) arms at 2.7 kpc from the nu- 
cleus values of ^ close to 80 km s _1 kpc -1 
would reduce shear effectively. He found ve- 
locity gradients close to 60 km s _1 kpc -1 for 
the NE arm and about 100 km s _1 kpc -1 
for the SW arm. For both arms, shear time 
scales were thus increased by about a factor of 
5 compared to the differential rotation case, 
with the shear becoming slightly reverse in 
the south. 

However, for the same regions, our con- 
clusions are different. In general, the gradi- 
ents that we measure are at least a factor 
of two higher than what Rand (1993b) finds 
for the southern arm crossing, implying that 
a condition of strong reverse shear will ex- 
ist. To judge from the above relation, this 
reverse shear will be even greater than shear 
in the pure differential rotation case. Shear 
timescales for selected GMAs and intercloud 
regions are listed in table 3, together with 
the corresponding Oort constants. In the re- 
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gions of extremely steep gradients found in 
between GMAs and through A10, this effect 
becomes so pronounced that it results in sig- 
nificantly shorter shear timescales than from a 
flat rotation. Thus, the conclusion that shear 
is generally reduced by streaming motions in 
spiral arms must be reconsidered. Of the three 
GMAs studied on the SW arm, only A9 shows 
reduced shear. Thus, conditions close to solid 
body rotation do occur, but apparently not 
throughout the arms. Future statistical stud- 
ies of cloud stability on the arms will have to 
take the varying degree of shear into account 
instead of generally assuming that shear is 
low. For more typical spirals where the aver- 
age streaming motions are smaller, the shear 
may still be generally reduced in the arms. 
Since the situation is likely to be complex, 
however, detailed high resolution studies of 
the velocity fields of the arms in a number of 
galaxies are necessery to decide this question. 

Elmegreen (1988) points out that shear 
may not be such an important factor in form- 
ing complexes when the gas is highly dissipa- 
tive, since the incident kinetic energy will be 
removed by processes such as cloud-cloud col- 
lisions. In this context, the variation in steep- 
ness of the velocity gradient along the arm 
may be interpreted as the GMAs having shed 
some of the strong reverse shear induced by 
the density wave, which still persists between 
the clouds. A10 may still be in the process of 
such a dissipation. It is quite suggestive that 
the region of reverse shear exists upstream of 
A9, suggesting the possibility that A9 once 
passed through this high-shear region on its 
way downstream. 

4-2.2. Interarm streaming 

The existence of the SI and S2 spurs makes 
it possible to study the velocity structure of 



gas in the interarm region. Figure 6c shows 
a position- velocity cut (pV-diagram) through 
the SI spur. The agreement with the RS pre- 
diction is good; the tangential velocity de- 
creases just before the gas enters the arm 
and then increases steeply in the arm. The 
expected peak in dispersion (Roberts 1993) 
right before the arm is also present in our 
data. We can measure the magnitude of in- 
terarm tangential streaming in these spur re- 
gions, and find that the gradient IS ~ 100 km 
s _1 kpc -1 , in the opposite direction to the arm 
gradient as qualitatively expected. This re- 
versed streaming is also expected at the out- 
side of the arm. 

There is a strong concentration of gas 15"- 
20" away from the M2 arm. This possible 
GMA seems to exist where the steep reversed 
tangential streaming begins to flatten out. 
Note that the dispersion in the gas also de- 
creases in this region. Perhaps the D clouds 
have formed in similar calm regions in spurs 
between the arms. 

Single dish studies suggest arm widths of 
20" - 30" instead of the 5" - 15" we find 
here. The structure along SI indicates that 
the velocity reversals occur 7" — 17" away 
from the arm center. The broader arm widths 
thus imply that the single dish surveys are 
picking up emission from the velocity rever- 
sal regions. The reversed streaming in these 
regions may be more disruptive to GMAs 
than even the high shear regions in the arm 
since the amount of dissipation will decrease 
with distance from the arm center. Tidal 
forces and flow expansion are other disrup- 
tive agents that may tear GMAs apart in the 
velocity reversal regions (Rand 1993b). Per- 
haps the CO emission here largely originates 
in diffuse unbound gas that the interferometer 
will filter out, but the single dish instruments 
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can detect. 

4.3. Radial streaming 

Figure 6d shows a pV-diagram along the 
minor axis which goes through the two main 
spiral arms at a radius of 25" to 35". Figure 6e 
shows a close up of the velocity discontinuity 
across the Ml arm. There is good agreement 
with the RS radial streaming model (dashed 
line in figure 6e), with the steep decrease in 
velocity, the density peak of the arm occuring 
at the velocity minima, and then the increase 
in velocity on the outer side of the arm. The 
magnitude of the shift is roughly 95 km s _1 . 
High resolution is needed to resolve the shifts 
and even at 2. "5 resolution, the shift is still 
somewhat unresolved, but an estimate is that 
it occurs on roughly 5" (230 pc), which means 
that the gradient is 430 km s -1 kpc -1 . Rand 
(1993b), however, finds only a 60 km s" 1 shift 
in the ionized gas (Ha) at the same arm cross- 
ing. The reason for this discrepancy is not 
clear, but if there is a displacement between 
the molecular and ionized gas, the Ha might 
miss the region of steepest gradient. Extinc- 
tion may also hamper the interpretation of 
the Ha emission. 

The radial streaming motions are not as 
pronounced on the eastern (M2 arm) minor 
axis. Although a drop in velocity can be de- 
tected on the inside of the arm, it is of smaller 
magnitude, 50 km s _1 , than on the west side, 
and the expected rise in velocity on the out- 
side of the arm is absent. The interpretation 
is also complicated by the multiple velocity 
structure of the gas inside the arm, at ra- 
dius 0-23". The figure reveals two emission 
regions, one peaking at 440 km s -1 , the other 
at 490 km s -1 , although they are at the same 
radius along the same cut through the mi- 
nor axis. Vogel et al. (1993) also note that, 



for their Fabry-Perot Ha data, the streaming 
for the M2 arm where it first crosses the mi- 
nor axis is weaker than what is measured for 
the Ml arm on the corresponding position. 
They attribute this to the more complex gas 
response (as indicated by dust lane morphol- 
ogy) of the M2 arm close to the oval distor- 
tion. 

Although streaming motions can be ob- 
served along the Ml arm all the way be- 
tween the minor and major axis arm cross- 
ings, very little streaming can be detected on 
the M2 arm until the major axis crossing in 
the north. The reason for this lack of stream- 
ing in the first quadrant of M2 is unclear. 
It can hardly be attributed to complex ve- 
locity structure which is suggested to explain 
the lack of streaming on the minor axis (see 
above). A possible explanation is that the M2 
arm in this region is only half as thick as the 
Ml arm and may be too narrow for a display 
of the full velocity range of the arm stream- 
ing. We would then expect to see streaming 
in M2 in another tracer than CO, sampling 
a different region of the arm. Vogel et al. 
(1993), however, note the lack of streaming 
in M2 in their Ha map as well. Perhaps the 
reason can be found in an assymmetry in the 
matter distribution. 

4-3.1. Galactic shocks 

The shape of the radial and tangential ve- 
locity components in figure 5 shows sharp de- 
celeration and troughs, signifying the pres- 
ence of a galactic shock. In the early models, 
self-gravity was not included and Lubow et al. 
(1986) suggested that introducing self-gravity 
would smooth density and velocity profiles, 
removing the shock structures. Roberts (1993) 
however reports that when including self-gravity 
in the simulations the abrupt velocity struc- 
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ture remains, even though the density profile 
becomes somewhat smoother. Our maps of- 
fer, for the first time, sufficient resolution to 
resolve the structure in the molecular stream- 
ing motions. As we have seen above, both 
the radial and tangential velocity components 
show steep gradients (see figues 6a, 6b, and 
6e) in accordance with predictions by the RS 
models. The region IC1 (in figure 6a), for 
example, clearly delineates the shock struc- 
ture formed on the southern arm crossing of 
the major axis. Our data thus support the 
presence of shocks in the arms of M51. The 
fact that the shock structure and magnitude 
varies along the arms indicates that local ef- 
fects enter to modify the effects of large-scale 
dynamics. 

Single dish surveys (Garcia-Burillo 1993ab; 
Kuno and Nakai 1997) report a more grad- 
ual change of the velocity vector across the 
arms. Kuno and Nakai (1997) find that the 
deceleration of the radial velocity component 
continues for about 1 kpc (20"). Because of 
this smoothness, they conclude that galactic 
shocks do not occur in the arms of M51. At 
their resolution (16"), however, the sharpness 
of the velocity profiles will be smeared, since 
the gradient occurs on a scale of 1" — 5" and 
we suggest insufficient resolution in the single 
dish studies is the main cause of the discrep- 
ancy. 

4.4. Ultraharmonic resonance 

Lubow (1993) and Bertin (1993) discuss 
the behaviour of the cold gas component at 
resonances and how one may use observations 
of gas in galaxies to detect resonance loca- 
tions. Due to its dissipative nature, the gas 
behaviour at the resonances is expected to be 
quite different from that of the stars, mask- 
ing the location of the resonance or, in some 



cases, revealing it through conspicious fea- 
tures. Shu, Milione and Roberts (1973) sug- 
gested that a second compression of the inter- 
stellar gas, associated with an ultraharmonic 
n=2 resonance, may generate secondary spi- 
ral features. Patsis et al. (1997) studied the 
response of a gaseous disk to a spiral pertur- 
bation, and their models show arm bifurca- 
tions and secondary arms between the main 
arms. They suggest that these bifurcations 
are the sign of the 4/1 ultraharmonic reso- 
nance and could serve as good resonance in- 
dicators. Artymowicz and Lubow (1992, AL) 
also predict the launching of a second wave in 
the gas at the location of the lowest order ul- 
traharmonic resonance. The resonance may 
result in the termination, or severe depres- 
sion, of the main arm (e.g. Contopoulos and 
Gr0sbol 1986; AL). Contopoulos and Gr0sbol 
also expect larger deviations from circular or- 
bits near the resonance which should result in 
increased radial streaming motions in the gas. 

The arm bifurcations in our figure 1 
symmetrically located on opposite sides of the 
galactic center at a radius of ~ 71" (3.3 kpc at 
D=9.6 Mpc). They are very similar to those 
in the models by Patsis and coworkers, which 
leads us to suggest that the Bl and B2 cloud 
chains may be faint, patchy, secondary spiral 
arms caused by the 4/1 resonance. According 
to Patsis et al. the actual resonance occurs 
somewhere between the arm bifurcation and 
the termination of the arm. Rand (1993b) 
places the 4/1 resonance in M51 close to 3.9 
kpc based on an intensity minima in the main 
CO arm. He also finds large streaming mo- 
tions (125 km s" 1 ) right before the minima 
at 85" radius. We find similar velocity shifts, 
130 km s _1 , in the gas at the location of the 
two bifurcations (see figure 3). 

If the rotation curve is known, the loca- 
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tion of one resonance will predict the posi- 
tion of the other resonances (e.g. Binney 
and Tremaine, 1987). There are, for exam- 
ple, several studies that have aimed at identi- 
fying the location of the corotation resonance 
(Tully 1974; EES; Vogel et al. (1993); Garci'a- 
Burillo 1993), see table 4, and we can use it to 
predict the expected location of the 4/1 reso- 
nance. We adopt R=6 kpc as the position of 
the corotation resonance, although the total 
range in table 4 is 5.8-7.4 kpc. For a Mestel 
disk (i.e. a disk in which surface density is 
inversely proportional to radius resulting in a 
flat rotation curve, e.g. Binney and Tremaine 
(1987)), having corotation at 6 kpc will place 
the 4/1 resonance roughly at 3.8 kpc. The 
estimated positions for the ILR lies between 
1.4-1.7 kpc (Tully 1974b; Garcia-Burillo et al. 
1993, see table 4), and for a Mesteldisk this 
places the 4/1 resonance between 3.0 and 3.9 
kpc. EES suggests that the Outer Lindblad 
Resonance (OLR) lies at R—8 kpc which, for 
the same model, places the 4/1 resonance at 
R=3 kpc. 

4 ■ 4 ■ 1 ■ Two spiral patterns 

According to Contopoulos and Gr0sbol (1986) 
the 4/1 resonance should represent the limit- 
ing radius for (strong) spiral arm structure, 
and the main arms should stop not far be- 
yond. AL agree that the 4/1 resonance may 
result in strong losses, but the amplitude fur- 
ther out at corotation could still be large. 

In M51, the spiral pattern clearly contin- 
ues beyond the suggested location of the 4/1 
resonance. The CO arm is indeed suppressed 
close to the resonance (e.g. Rand 1993b, see 
above) but commences again at larger ra- 
dius. It has been suggested that there are 
two pattern speeds in M51, with the outer 
arms belonging to a different pattern, driven 



by the companion galaxy NGC 5195 (Tully 
1974b; Elmegreen, Elmegreen, and Seiden 
1989 (EES); Vogel et al. 1993). EES, for ex- 
ample, propose 1. An inner spiral mode that 
has an Outer Lindblad Resonance (OLR) at 
the position of a prominent arm intensity gap 
at R—8 kpc (pattern speed 90 km s -1 ), 2. 
and an outer material arm that corotates with 
the companion (pattern speed 22 km s _1 ). 
They argue that the position of the ILR of 
the outer arms occurs at the same radius as 
the corotation resonance of the inner spiral, 
and they suggest that the outer spiral trig- 
gers a response in the vicinity of its ILR which 
stimulates corotation of an inner spiral wave- 
mode. The radius of the ILR of the outer spi- 
ral mode represents the inner limit to its ex- 
tent and therefore, the inner spiral structure 
extends beyond its 4/1 resonance, at least out 
to corotation. 

Toomre (1981) argues that the grand de- 
sign spiral structure of M51 is largely the re- 
sult of a strong, transient perturbation, and 
that the density waves excited by it are likely 
to differ significantly from those considered in 
the density wave theory of isolated systems. 
Specifically, the 4/1 resonance should not be 
excited in such a transient encounter. As we 
have discussed , however, the arm bifurca- 
tions, CO intensity minima and large stream- 
ing motions are expected signs of the 4/1 res- 
onance. EES also discuss optical evidence for 
the presence of a 4/1 resonance. Further stud- 
ies, theoretical and observational, are neces- 
sary to resolve this issue. 

4.5. Molecular structure in the central 
arcminute 
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4-5.1. Spiral structure 

The Inner Lindblad Resonance (ILR) oc- 
curs at the apparent termination of spiral 
structure at a radius of 25" — 30" (e.g. Tully 
1974b; this work). However, narrow tongue- 
like features in the central 25" of our ro- 
bustly weighted CO map (figure 7) hint at 
the possibility of an inner compact pattern 
in M51. The possibility of such structure 
was first raised by Zaritsky, Rix and Rieke 
(1993, ZRR), who from a K-band image resid- 
ual found evidence of spiral arms inside the 
ILR. They claim that these arms continue to 
wind through an additional 540° beyond the 
ILR, ending at about 10" from the nucleus. 
They further suggest the existence in the in- 
nermost 10" of a small bar, possibly an oval 
bulge. It is surprising to find spiral arms in- 
side the ILR since it is generally claimed that 
a density wave will not propagate there (e.g. 
Lubow 1993). By analogy with the outer ma- 
terial arms, however, the structure within the 
ILR may not be part of the same pattern as 
the spiral structure beyond a radius of 25": 
arms, for example, may be driven by a central 
bar, or by the spiral pattern further out. It is 
also possible that the assumed position of the 
ILR is wrong and should be placed at much 
shorter radius. This is a possibility since it 
is difficult to measure the rotation curve in 
the very center of the galaxy. We have at- 
tempted to compare the ZRR K-band image 
of the inner spiral structure with our robustly 
weighted CO map to see if the spiral arms in 
the ZRR map (dashed lines in figure 7) can be 
identified with the features in the molecular 
gas. Some of the narrow CO features in the 
inner arcminute are found to agree in posi- 
tion with the K-band arms, but some do not. 
The width of these tentative molecular arms 
is about 5" and the structure is patchy, not 



tracing the arms fully, as is also true for the 
K-band arms. 

4-5.2. The central bar 

The presence of a bar in the inner region 
of M51 has been discussed by several authors 
(e.g. Tully 1974b; Smith et al. 1990; Pierce 
1986; ZRR; Garcia-Burillo et al. 1993b; Kohno 
et al. 1997). Tully (1974) identifies the 
ILR with the observed terminal radius of the 
(main) spiral structure at 25" and from de- 
viations of circular motions in his Ha data 
suggests the existence of an oval dispersion 
ring inside the ILR with PA 130°, causing 
gas streaming on elliptical orbits, as found by 
Tully (1974b). There seems to be a general 
consensus that the PA of the bar, if it exists, 
is approximately 130°, but the length of the 
bar is fairly uncertain: between 10" and 20" 
(e.g. Pierce 1986; ZRR). 

Our data show no continuous gas structure 
that fills the central bar. This may be because 
the bar is filled with diffuse faint molecular 
gas that the interferometer will filter out, or 
simply because the bar is deficient in molec- 
ular gas, or it may mean there is no central 
bar at all. The robustly weighted map does 
show some faint emission connecting the cen- 
tral gas structure to the emission region to 
the northwest (see fig lb) but this structure 
is only a few arcseconds long; discussion of it 
will be left to paper II. 

In order to look for evidence of non-circular 
motions in the center, we used the AIPS GAL 
task to subtract a rotation curve from the nat- 
urally weighted velocity field. Initial values 
were adopted from Tully (1974a): PA=170°, 
i = 20°, v sys =464 km s _1 . We adopted a rota- 
tion curve of the form -2- = l-e~ ln{W0 T^ ) . 
The resulting central residual map is shown in 
figure 8. In the inner 30" there are two main 



18 



deviations from the general rotation curve, a 
blueshifted region to the east and a redshifted 
region to the west. The size of the region of 
deviating velocities is roughly 20" x 10" and 
the PA of the joined structures is about 125°. 
The magnitude of these deviations is 20-30 
km s" 1 and they are consistent with gas on 
elliptical orbits in a bar. The region of large 
residual velocities agrees roughly with the end 
of the K-band spiral arms of ZRR - - and 
with the beginning of the bar they suggest 
is present in the inner 10" — 20". 

4-5.3. The nature of structure C 

The structure we have named the C fea- 
ture (figure la) is located at the inner end 
of the main NE (M2) spiral arm, at what we 
agree is the radius of the ILR of M51. The fea- 
ture also lies at the eastern end of the possible 
central bar, although we note that the corre- 
sponding feature at the western end is not so 
spectacular. Tully points out that the regions 
where the spiral arms meet the material on el- 
liptical orbits (the "dispersion ring") usually 
show evidence of abrupt velocity discontinu- 
ities, pronounced dust lanes and strong radio 
continuum emission. In our data there are 
multiple spectral peaks at several positions in 
this region. As we can see in the dispersion 
map (figure 4), the very inner region of the 
C feature has quite narrow lines (low disper- 
sion) while the more complex velocity struc- 
ture occurs around the "edges" of the feature 
(as indicated by red regions in figure 4). In 
the lower (southern) region there are multi- 
ple spectral featues, where the less prominent 
feature is at velocities 20-50 km s _1 lower 
than the main feature (figure 9a). It is pos- 
sible that the low-velocity-feature (which is 
about 3 times fainter than the main one) orig- 
inates in gas on bar orbits while the brighter 



emission comes from gas on circular orbits. 
The upper (northern) "red" region also con- 
sists of double spectral components, but here 
the high velocity feature is the fainter one (fig- 
ure 9b). There are also some positions where 
the high dispersion is caused by broad lines, 
> 50 km s -1 , rather than multiple spectra. 
We speculate that the disturbances in feature 
C are the result of high velocity collisions be- 
tween molecular clouds in different orbits. If 
cloud-cloud collisions occur anywhere, they 
are most plausible in such a region, where cen- 
tral spiral arms, main spiral arms, and gas on 
elliptical orbits converge. 

5. Conclusions 

We have undertaken a new sensitive, high res- 
olution CO study of the inner 2/5 of M51. 
The main results are: 

The emission in the two major arms 
(NE and SW) originates mostly from 
supermassive cloud complexes (Giant 
Molecular Associations, GMAs) which 
are for the first time resolved both along 
and perpendicular to the arms. The 
brightest (T B (CO) = 8 K) GMA on the 
SW arm has a FWHM size of 11" by 
6" (512 by 280 pc) and sits on a ridge 
of CO emission 20" x 7" in size. Along 
the SW arm we count 16 GMAs, with 
sizes varying from 6" x 3.5" to 14" x 8". 
The overall morphology of the CO emis- 
sion is symmetric in reflection about the 
nucleus with major complexes occuring 
opposite each other in the two major 
arms. The CO line flux in the inner 
area of the SW arm is approximately 
twice as bright as that from the analo- 
gous location in the NE arm. The typ- 
ical one- dimensional velocity dispersion 
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is 7-15 km s _1 on the arms and the arm 
width varies between 5"- 15". The SW is 
broader than the NE arm by ~50%. 

The average aperture synthesis arm-to- 
interarm CO flux ratio is about 12. We 
detect five times more flux in the arms 
than in the interarm region. The to- 
tal flux detected by the interferometer is 
2.7 x 10 3 Jy km s" 1 , which is about 50% 
of the Onsala 20m single dish flux in the 
same region. For a standard Galactic 
CO mass calibration factor (see table 1) 
the molecular mass seen by the interfer- 
ometer is 1.9 x 10 9 M , 15% of which 
lies in the interarm regions. 

Streaming motions due to density waves 
are evident along the minor and major 
axis. The magnitude of the streaming 
motions is generally very high, 60-150 
km s" 1 . Streaming, however, is ab- 
sent on the NE arm beginning at the 
eastern minor axis until the arm crosses 
the major axis. Our maps offer, for 
the first time, sufficient resolution to de- 
tect structure in the molecular stream- 
ing motions. Both the radial and tan- 
gential velocity components show steep 
gradients in qualitative accordance with 
predictions by the density wave mod- 
els of Roberts and Stewart (1987). The 
region IC1 (in figure 6a), for example, 
clearly delineates the shock structure 
formed on the southern arm crossing of 
the major axis. Our data thus support 
the presence of shocks in the arms of 
M51. 

Streaming motions (at least the tangen- 
tial ones) vary in structure along the 
spiral arms. The velocity gradients be- 
come generally steeper between GMAs 



than within them; gradients in places 
reach 1000 km s _1 kpc -1 over a small 
region of 100-200 pc. Velocity gradients 
in GMAs are quite variable from only 
60 km s" 1 kpc~ x in one GMA to sev- 
eral hundred km s _1 kpc -1 in another. 
Therefore, the notion that shear is gen- 
erally reduced by streaming motions in 
spiral arms will have to be modified. Of 
the three GMAs studied on the SW arm, 
only one enjoys a shear-reduced envi- 
ronment. Conditions close to solid body 
rotation do occur, but apparently not 
everywhere in the arms. Future statisti- 
cal studies of cloud stability on the arms 
will have to take the varying degree of 
shear into account instead of a general 
assumption that shear is generally low. 

There is evidence of bifurcations in the 
molecular spiral arms. Emerging from 
these bifurcations are two chains of molec- 
ular clouds, Bl and B2. These features 
may be evidence of a secondary com- 
pression in the density wave caused by 
the ultraharmonic 4/1 resonance. We 
suggest that the arm bifurcations occur 
near the 4/1 resonance of the spiral pat- 
tern, which is then placed between 3.2 
and 3.8 kpc from the galactic center. 
This location of the resonance is consis- 
tent with results by Rand (1993b). The 
4/1 resonance is well inside the M51 spi- 
ral pattern and this supports the sug- 
gestion that M51 has at least two spiral 
systems. 

Other interarm features found are a 
short chain of low- dispersion clouds mid- 
way between the main and secondary 
arms and several molecular spurs pro- 
truding from the main spiral arms, in 
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particular on the western side of the 
map. Tangential interarm streaming 
motions are detected along two of the 
most prominent spurs. The agreement 
with the RS prediction is good; the tan- 
gential velocity decreases just before the 
gas enters the arm and then increases 
steeply in the arm. The gradient of the 
interarm tangential streaming is ~ 100 
km s _1 kpc -1 . 

There is evidence in our highest resolu- 
tion map (with robust weighting) for a 
central molecular spiral inside the claimed 
location of the ILR; it is consistent with 
the K-band spiral found by Zaritsky, 
Rix and Rieke (1993). 

We find evidence of non-circular mo- 
tions in the inner 20" of M51. The mag- 
nitude of these deviations is 20-30 km 
s _1 and they are consistent with gas on 
elliptical orbits in a bar. Our data, how- 
ever, show no continuous gas structure 
that fills the central bar. 

We find an unusual structure, an expan- 
sion in the NE molecular arm at 25" ra- 
dius (1.2 kpc) southeast of the center. 
This broadening occurs right after the 
end of the NE arm at the Inner Lindblad 
Resonance (ILR). Multiple-peak spec- 
tra, velocity twists and structures with 
apparently high velocity dispersion are 
associated with this feature. We spec- 
ulate that the disturbances in feature 
C are the result of high velocity colli- 
sions between molecular clouds in dif- 
ferent orbits. 
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Fig. 1. — it a) CO 1-0 integrated inten- 
sity map, naturally weighted with synthesized 
beam 3"95 x 3"27 and BPA=-42.8°. Contour 
levels are 1.1 Jy beam -1 km s _1 x (1.0, 2.5, 
4.0, 17.5, 19.0). The peak flux is 29.04 
Jy beam" 1 km s" 1 at a = 13 /l 27 m 43 s .715 ; 
5 = +47°26'52".25 in the Ml arm. The low- 
est contour is at the 3a level. The total in- 
tegrated flux in the map is 2.7 x 10 3 Jy km 
s" 1 . The solid lines mark the secondary arms 
Bl and B2 and the dashed lines mark the 
string of interarm D clouds (see text). Ar- 
rows indicate other features discussed in the 
text. The outer map cutoff is indicated with 
a dashed line. This cutoff is approximately 
5" inside the maps outer primary beam half 
power points, b) CO 1-0 integrated inten- 
sity map robustly weighted with synthesized 
beam 2."88x2."ll and BPA=-80.7°. Contour 
levels are 0.4 Jy beam" 1 km s" 1 x (1.0, 4.0, 
8.0, 32). The peak flux is 18.64 Jy beam" 1 
km s" 1 at the same position as for the nat- 
urally weighted map. The lowest contour is 
at the 3a level. The total integrated flux in 
the map is 2.0 x 10 3 Jy km s" 1 . The outer 
map cutoff is indicated with a dashed line, c) 
Position- Velocity cut through the Ml and Bl 
features, at PA 242°. Zero velocity is at 472 
km s" 1 . 

Fig. 2. — a) The CO radial brightness distri- 
bution for the naturally weighted map which 
shows the brightness for each radial bin out 
to a radius of 95" (corrected for inclination). 
For the conversion factor used in this paper, 
Jy beam -1 km s" 1 translates to gas surface 
density, b) Intensity cut along the Ml and 
M2 arms. The resolution has been smoothed 
to 10". 

Fig. 3. — The CO 1-0 velocity field on the 
naturally weighted map. The scale range 



from 350 km s 1 (blue) to 550 (red) km s 1 , 
the contour spacing is 12.5 km s" 1 . 

Fig. 4. — The CO 1-0 one-dimensional dis- 
persion map. The scale ranges from 4 (blue) 
to 22 (red) km s" 1 . 

Fig. 5. — The cloud density (upper panel) 
and velocity componets (tangential (middle 
panel) and radial (lower panel)) plotted vs. 
spiral phase, from RS dissipative simulations, 
case M. For a tightly wound spiral, a cut 
through spiral phase is roughly equivalent to 
a cut in radius. This particular figure is at a 
time step 1170 Myr. 

Fig. 6.— a) Velocity field (2. "5 resolu- 
tion) overlayed on the integrated intensity 
map (grayscale) in the region where arm Ml 
crosses the major axis. Cloud complexes 
A8,A9,A10 are indicated in the figure. The 
contour spacing is 5 km s" 1 , the grayscale 
ranges from to 12 Jy beam" 1 km s" 1 . b) 
Velocity field (2. "5 resolution) overlayed on 
the integrated intensity map (grayscale) in 
the region where arm M2 crosses the major 
axis. Cloud complexes Al, A2 are indicated 
in the figure. The contour spacing is 5 km 
s" 1 , the grayscale ranges from to 8 Jy 
beam" 1 km s" 1 . c) Position- Velocity cut 
through the Al GMA and the SI spur fea- 
ture, at PA 90°. The velocity decreases (on 
the blueshifted side of the galaxy this trans- 
lates to a redshift) right before the gas enters 
the arm and then increases (appears as a de- 
crease on the blueshifted side) steeply in the 
arm. The expected (Roberts 1993) peak in 
dispersion right before the arm is also present 

d) Position- Velocity cut along the minor axis. 

e) Blow-up of the previous figure showing only 
the cut through the western Ml arm. In- 
dicated as a dashed line is the Roberts and 
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Stewart (1987) model of radial streaming mo- 
tion. 

Fig. 7. — Blow-up of the central region of the 
robustly weighted moment map, figure lb. Ri- 
band spiral arms are indicated with dashed 
lines. 

Fig. 8. — Inner region of the residual velocity 
map after a rotation curve has been removed. 
Solid lines indicate redshifted velocities with 
respect to the model = l-e"'^ 100 *^). 

^ V V max ' 

Contour spacing is 10 km s" 1 . Grayscale 
goes from 10 to 40 km s" 1 residual veloci- 
ties. Regions of blueshifted deviations from 
the model are indicated with patterns: verti- 
cal lines indicate velocities t — 30 km s _1 ; 
checks indicate -30 — -20 km s _1 ; horizontal 
lines indicate -20 — -10 km s" 1 . White in- 
dicates regions of no deviation from circular 
velocities. The greatest deviations from circu- 
lar rotation are found in two 10" x 20" regions 
west (redshifted) and east (blueshifted) of the 
center; the small spurs (or, rather fractions of 
spurs) protruding from the western arm also 
shows deviations from circular velocities by 
up to 40 km s™ 1 ; there is also a redshifted 
region just north of the blueshifted region to 
the east. The nuclear disk is redshifted with 
respect to the systemic velocity. 

Fig. 9. — a) Sample spectrum from the 
southern edge of the C feature at position 
a: 13 h 27 m 47 s .608, 8: 47°26'50"25. b) Sam- 
ple spectrum from the northern edge of the 
C feature at position a: 13 h 27 m 47 s .806, 5: 
47°27'04"25. 
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Table 1 
M51: Adopted Properties 



Parameter Value 



Center position (1950.0) a a: 13 h 27 m 46 s .327 

5: 47°27'10"25 

Morphological type b Sbc 

Systemic Velocity (hel) c 464 km s _1 

Distance d 9.6 Mpc 

Position angle c 170° 

Inclination 20° 

Linear resolution on galaxy (robust) 134 x 98 pc 

Linear resolution on galaxy (natural) 184 x 152 pc 

Adopted conversion factor JV(H 2 )/I(CO)=2.3 x 10 20 cnr 2 (K km s" 1 y 1 



a Position of radio continuum peak (Ford et al. 1985) 
b Sandage and Tammann 1981 
c Tully 1974a 

d Sandage and Tammann 1975 

C A standard Galactic 12 CO luminosity to H 2 mass ratio (Strong et al. 1988). For 
observed CO fluxes, this translates to M(H 2 ) = 0.84 x lO^SAvD 2 M (D is the distance 
in Mpc; SAv is the integrated 12 CO 1-0 line flux in Jy km s _1 ) 
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Table 2 
M51: Map Pointing Centers 



Pointing 
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a Map center position 
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Table 3 

M51: Tangential streaming and shear 



Region ^ A(R) h t c shear 

(km s _1 kpc -1 ) (km s _1 kpc -1 ) (10 7 yr) 

GMAs: 

A2 250 -85 1.1 

Al 115 -18 5.3 

A8 250 -85 1.1 

A9 65 8 12 

A10 700 -310 0.3 

Inter-cloud regions: 

IC1 800 - 1300 -360 — -610 0.3 — 0.2 

IC2 600 -260 0.4 

IC3 400 -160 0.6 



a Velocity gradient in the tangential streaming motions 

b Oort constant at radius R, A(R) = \{j[ — ^) where v c is the circular 
velocity. At R=2.7 kpc ^ is 80 km s -1 kpc -1 . 

Tor pure differential rotation, t s hear is 2.4 x 10 7 yr at R—2.7 kpc. 
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Table 4 
M51: Position of resonances 



Reference 13 ILR corotation Comments 

Garcia-Burillo 1 1.7 7.4 assumes single pattern 

Tully 2 1.4 5.7-7.1 two spiral patterns 

Elmegreen 3 • • • 6.0 two spiral pattern 

Vogel 4 • • • 6.0 two spiral patterns 

a All data has been corrected to a distance .0=9.6 Mpc 

b (l) Garcia-Burillo et al. (1993b); (2)Tully 1974b; (3)Elmegreen, Elmegreen and Seiden (1989); (4) 
Vogel et al. (1993) 
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Radial CO brightness distribution 
(surface brightness distribution) 
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